Resonant Manipulation of d-wave Interaction of Cold Atoms with Two Lasers and a Magnetic Field 
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We present a theory for manipulation of d-wave interaction of cold atoms with two lasers strongly driving 
two photoassociative transitions. The theory predicts the occurrence of a coherence between two excited ro- 
vibrational bound states due to the photoassociative dipole-couplings of ground-state d-wave scattering state 
to the bound states. We show that this excited-state coherence significantly influences atom-atom interaction. 
In particular, this leads to the enhancement of d-wave elastic scattering and to the suppression of inelastic 
scattering. In the presence of an s-wave magnetic Feshbach resonance, the two lasers can couple the s-wave 
resonance with the d-wave scattering state leading to the further enhancement in d-wave scattering at relatively 
low energy. Our numerical calculations based on realistic parameters show that d-wave manipulation would be 
most effective in case of atoms having excited diatomic states with narrow natural linewidth. We estimate that 
at 100 pK the inelastic scattering rate in Yb can be reduced to 20 s _1 while the elastic scattering rate can be 
two orders of magnitude larger. 

PACS numbers: 67.10.Db,03.65.Nk,37.10.De,74.90.+n 



I. INTRODUCTION 

Cold atoms have now become a testing ground for inter- 
acting many-particle physics, in particular for models of con- 
densed matter systems. They offer a unique opportunity to 
study many -body physics lU with a controllable two-body in- 
teraction. By making use of the tunability of the s-wave scat- 
tering length with an external magnetic field, s-wave physics 
of interacting ultracold atoms has been extensively studied in 
recent times. To go beyond the s-wave physics, it is essential 
to devise methods for manipulation of higher partial-wave in- 
teractions of cold atoms. Higher partial-wave interactions are 
important for some exotic phases of matter such as the p-wave 
anisotropic superfluidity of liquid He-3 and the d-wave su- 
perconductivity of cuprate materials Q2jQ . The ability to study 
phenomena related to d-wave Cooper pairing with fermionic 
atoms will help us to gain new insight into the open problem of 
high temperature superconductivity. To explore d-wave pair- 
ing of fermions with two spin components (like electrons) at 
various interaction regimes, it is of primary interest to devise 
an effective method of manipulating the d-wave interaction in 
two-component fermionic atoms. 

There are mainly two methods of controlling atom-atom 
interaction at ultracold temperatures. The most widely used 
method involves the magnetic Feshbach resonance (MFR) 
JH 0]. The other method using the optical Feshbach reso- 
nance (OFR) JH 0] is a relatively recent one. Over more 
than a decade, the MFR has been used to study the s-wave 
physics of interacting atomic gases. In fact, the MFR has 
become a standard tool for studying Fermi superfluidity [jvj] 
and strongly interacting Fermi gases 13] with tunable s-wave 
scattering length. The OFR has been originally developed 
as a method for altering the s-wave scattering length as in 
the MFR. This method relies on off -resonant free-bound pho- 
toassociative transitions. The use of strong photoassociative 
coupling has been proposed as an optical method of manip- 
ulating higher partial-wave interactions [9]. Of late, p-wave 
OFR has been experimentally reported in fermionic 171 Yb 
atoms by Yamazaki et al. IU0I1 based on a suggestion by 



Goyal et al. Hill to use purely long-range states lfl2l [T3ll 
of photoassociated excited molecules. Earlier, magnetic field 
induced p-wave Feshbach resonances have been observed in 
spin-polarized 40 K ifTill and 6 Li lfl5ll Fermi gases. Not only p- 
wave OFR, but alo generating other higher partial-wave inter- 
actions in ultracold atoms under various physical conditions 
has attracted a great deal of research interest lUrjl - COl in recent 
times. The generation of d-wave interaction in ultracold atoms 
by strong confinement has been discussed in Refs lfl8l |20ll . 
However, an efficient manipulation of the d-wave atom-atom 
interaction for the purpose of studying d-wave many-body 
physics with cold atoms is yet to be demonstrated. Both p- and 
d-wave scattering amplitudes of cold atoms are in general ex- 
tremely small according to Wigner threshold laws. This poses 
a major challenge in manipulating p- and d-wave interactions 
with an optical or any other field. The method we propose 
here is applicable for both p- and d-wave manipulation, but in 
this work we specifically concentrate on d-wave only. 

Here we present a theory for an all optical or magneto- 
optical manipulation of d-wave interaction in three dimen- 
sions. Consider two colliding ground-state cold atoms irradi- 
ated by two photoassociation (PA) lasers tuned near resonance 
with two excited ro-vibrational bound (molecular) states 1 and 
2 characterized by two different rotational quantum numbers 
J 1 and J r 2, respectively. However, these states may have the 
same or different vibrational quantum numbers. These ro- 
vibrational states must be accessible by PA transitions from 
the d-wave scattering state of the two ground-state atoms. The 
method we propose here relies on the quantum interference 
between the two PA transition pathways leading to the occur- 
rence of a coherence between the two excited states. We show 
that this coherence causes suppression of inelastic scattering 
and enhancement of d-wave elastic scattering. Note that the 
inelastic scattering considered here is due to absorption from 
the PA lasers and loss due to spontaneous emission. Because 
of spontaneous emission from excited states, one-color OFR is 
bound to suffer from inelastic loss as discussed by Bohn and 
Julienne I2H1 . However, here we show that, by establishing 
a coherence between two excited ro-vibrational levels by two 
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FIG. 1: Plots of s- and d-wave elastic scattering cross sections ao 
(dashed-dotted curves) and 02 (solid curves), respectively, and total 
inelastic scattering cross section a ine i (dashed curve ) for Li atoms 
as a function of dimensionless scaled intensity (f/2 /7sp) of laser 2 
for different values of the intensity of laser 1 and magnetic field: (a) 
gi - 0.5-y 3p and (B - Bo)/|A| = -0.24, ( where A is the width 
of the MFR) (b) pi = 0.57 sp and B > B , (c) gi = 107 sp and 
(B - Bo)/ A = -0.24, and (d) g x = 107 sp and B > B . 7 sp is 
assumed to be 2 MHz and energy E = 20 fjK. 

lasers, it is possible to suppress this inelastic loss significantly. 
Previous studies ||22T - I24ll on Fano-Feshbach resonances in PA 
l22tl in the presence of a MFR j25l - l2~7ll primarily focused on 
controlling s-wave scattering states, suppression of inelastic 
scattering rf23tl and control over power broadening [l24ll . The 
present work is basically concerned with optical manipulation 
of the anisotropic atom-atom interaction. In case of magnetic 
atoms, a magnetic field may be applied along the z axis to in- 
duce an s-wave MFR. For two-component fermionic atoms or 
spin-polarized bosonic atoms, it is possible to optically couple 
the s-wave MFR to the rf-wave scattering state by exploiting 
appropriate rotational selection rules for the two photoasso- 
ciative transitions as analyzed below. The results of our nu- 
merical calculations based on realistic parameters reveal that 
the optical manipulation of the d-wave interaction is most ef- 
fective in atoms having narrow natural linewidths such as Yb. 

H. ANISOTROPIC DRESSED CONTINUUM 

Here we present a dressed state description of the ground- 
state continuum interacting with a magnetic field and a pair of 
laser fields. 



A. The model and its analytical solution 

For simplicity, let us consider only two ground-state scat- 
tering channels | g\) and | g 2 ) with | g\) being open (contin- 
uum) and I 172} closed. To keep our discussions most general 



at the outset, we here do not specify the compositions of the 
channels, but we will give details of the channels in a specific 
numerical illustration later. The closed channel is assumed 
to support an s-wave bound state | x) whose binding energy 
lies close to the threshold of the | <?i) channel. These two 
ground-state channel thresholds and the closed-channel bound 
state may be tunable with an external magnetic field leading 
to an s-wave Feshbach resonance. The ground-state channels 
are asymptotically diagonal. However, at intermediate sepa- 
rations there is an interchannel coupling V(r) which depends 
on the interatomic separation r. 

The energy spacings between the excited bound states | bi) 
and I 62) characterized essentially by the two different rota- 
tional quantum numbers J\ and J2 are assumed to be large 
enough so that the two lasers L\ and L 2 can drive only these 
two bound levels, respectively. If these two bound states be- 
long to the same vibrational level v, then the rotational energy 
spacing between them must be large enough. Furthermore, 
these bound levels are assumed to be far below the respec- 
tive dissociation threshold so that the transition probability 
at the single-atom level is negligible. The excited potentials 
asymptotically correlate to two separated atoms with one be- 
ing in the electronic ground 5 state and another in the excited 
P state. If the two atoms are homonuclear, then the excited 
state potentials asymptotically go as C^/r 3 where C3 repre- 
sents the resonant dipole-dipole interaction coefficient. In the 
case of heteronuclear atoms, the excited-state potential goes 
as ~ — 1/r 6 at large separations. 

In the rotating wave approximation, the Hamiltonian of 
these three bound states interacting with the ground contin- 
uum can be expressed as H = Hq + Hi where 

H = ^2{E n -fuv Ln ) I b n )(b n I 

n 

+ E X \ X )(X I ® I 22X52 + E/ E ' dE ' 

x I E' lmi,) harcharc (E' tmi | <g> | gi)(gi | (1) 

is the free part of the Hamiltonian with E n representing the 
binding energy of n-th excited molecular state | b n ), E x be- 
ing the energy of the cosed-channel bound state | \) an d I 
E', imi )baro being the bare (unperturbed) partial-wave (Imi) 
scattering state with energy E '. Note that all the energies are 
measured from the open-channel threshold unless stated oth- 
erwise. The interaction part of the Hamiltonian is 

Hi = EE/ dE '^TSiSE') I b n ){E'lm t \ ( 9l \ 

n,M„ lm e 

+ J dE'V x0 {E') I X ) barc<£'00 I ® I g 2 )( 9l I 

+ n "x I b n)(X I (92 I +C.C. (2) 

n,M n 

Here " (E) is the dipole matrix element of transition | 
b n ) EJm e } halc I gi), V x q(E) is the coupling between 
the quasibound (x) state and the s-wave (£ — 0) scattering 
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state of the bare continuum and fl nx is the Rabi frequency 
between the nth excited bound state and the \ state. 

Let | E, k) represent an eigen ket that satisfies the time- 
independent Schrodinger equation H \ E, k) = E \ E, k), 
where k denotes a unit vector along the incident relative mo- 
mentum of the two atoms. These eigen kets are energy nor- 
malized in a manner such that (E k \ E,k) = 6{E — 
E')5(k — k 1 ). To solve the Schrodinger equation, we expand 
| E, k) in the following form 

I E,k) = Y, Y AnE I 6 " 

(J n ,M n )) + B E I x) .92) 

n=l,2 M n 

E E / dE 'C l E '%Yi'ra t Xk) I E',em,) haie I .91) 



l,m£ £' ,rrig 



(3) 



where Be, A„e and C E 7 l e t ^ ne are the expansion coefficients 
to be derived. The coefficients Be and A u e can be ex- 
panded as B E = Efm ( ,4 m ' T fm ( ^) and A nE = 

Et> me , A nE e 'Yi' mt M Physically, <>' implies the prob- 
ability amplitude for the excitation of the bound state | b n ) 
when the incident partial-wave of relative motion of the two 
atoms is I' and its projection along the space-fixed z axis is 
my. Similarly, B E t denotes the probability amplitude for 
the occupation of the quasibound state | x) f° r the (i'my) 
incident partial wave. The analytical expressions of the coef- 
ficients A e nE ' ' , B E me ' and C E 7 l e e ^ ne are derived in Appendix 
A. They are expressed in terms of the basic coupling param- 
eters of the model, namely the laser couplings Aj" 1 ^ and 
magnetic coupling V x o. Explicitly, 
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$n + l*fi(Sfn + 7n)/2 ' 



(4) 



The first term &/ n o in the numerator of Eq. (|4]i is defined 
by Eq. ( |A25t in Appendix A. This term arises from inci- 
dent partial wave £' = (s wave) only. It consists of two 
terms. The first term /3„ e Aj° M is the amplitude for transi- 
tion form the s-wave scattering state to the nth bound state 
in the presence of a MFR and laser L n . The effect of the 
MFR is given by the dimensionless parameter f3 ne = qn S^ e . 
The parameter q„ f is the Fano-Feshbach asymmetry parame- 
ter 12411 . e is a dimensionless energy parameter defined by Eq. 
(lA19t and related to the MFR phase shift r) ies by the relation 
cot r] res = —€. In the limit e — > ±00 or equivalently r) res — > 
or the coupling V x o — > we have (3 n —> 1, The second term 
7 /C n „'/3„' £ A , M on the right-hand side (RHS) of Eq. 

(1A25I > describes the combined effects of the other (n'th) bound 
state, two lasers and the MFR on the transition amplitude from 
the s-wave scattering state to the n-th n') bound state. The 
quantity K nn , = + K™„ where K%£, defined by Eq. 

(IA23b in Appendix A is an effective coupling between the 
two excited bound states induced by the two lasers in the pres- 
ence of the MFR provided both the bound states are accessible 
by PA transitions from the s-wave scattering state. IC™^, de- 
pends on the two Fano-Feshbach asymmetry parameters q\ / 



and 92/ and therefore on the optical dipole transitions from 
the closed channel s-wave quasibound state to the two excited 
states. /C^/ given by Eq. (IA17b in Appendix A represents 
an effective coupling between the two excited states induced 
by the two lasers through continuum-bound PA transitions. 
Note that fc^i depends on the properties of the open-channel 
bare continuum and does not depend on closed-channel qua- 
sibound state or the bound-bound transition. All possible in- 
cident partial waves which are allowed by selection rules for 
both the PA transitions contribute to /C2^, . The second term 
&nt> in the numerator of Eq. (|4]i is given by Eq. (1A261 in 
Appendix-A. This term describes transition amplitude from a 
nonzero incident partial wave (£' > 0) to the nth bound state 
in the presence of both the lasers and the MFR. In the denom- 
inator of Eq. ©, S n = E - M n - Ef m - g£* - S ml ,, 
where E~ is the light shift of n-t\\ bound state due to 
L n laser only. given by Eq. JA33t in Appendix A 

is a shift that strongly depends on the MFR, in the limit 
e —> ±00 this shift vanishes. S nn i defined by Eq. (IA341 > in 
the Appendix-A depends on both lasers through the coupling 
K nn i . The total linewidth for the transition to the nth bound 
state is 5f„ + j n where j n is the spontaneous linewidth and 
% = r "oo+X^>o,m £ T nem t +Tnn' depends on the magnetic 



and laser fields. Here T™q 



| 2 r„oo which, in the limit 



e — > ±00 reduces to the s-wave partial stimulated linewidth 
r„oo- Thus in the absence of the MFR, the first two terms in 
'Sn can be added to yield r„ = J^e m ^nimt which is the 
stimulated linewidth of the nth bound state due to laser L n . 
T nn r is defined by Eq. dA36b in Appendix A and explicitly 
depends on both lasers through K, nn i . 

Finally, the remaining expansion coefficients of Eq. (|3), the 
detailed derivation of which is given in the Appendix-A, can 
be expressed as 



B 



(e + *)r, 



V xQ {E)8 V0 + J2 (inf-i) 



x nA^V x0 (E)A^' 



Ce™^(E) = S(E - E')8 u .8 mtmt , + 



(5) 



Ff™/(E,E>) 



E — E' 



(6) 



where 



F e e 'Z l '(E,E') 



= B 



E 



'V x0 (E')S e , 



/ , ^nE ly lm t 



(E>). 



(7) 



Equations ©, ©, © and (0 clearly show that all the ex- 
pansion coefficients of Eq. © can be evaluated in terms of 
the known input parameters of our model. The effects of the 
two lasers and the magnetic field on atom-atom scattering are 
given by Eq. (0 as discussed in the next section. 



4 



B. Elastic and inelastic scattering 

Making use of the analytical solution as outlined above, 
we now derive scattering T and S matrices. The asymptotic 
form of the dressed continuum wave function (r | E, k) ~ 
\E , £;(fc,r), where r denotes relative coordinate, is given by 

$ E (k, r oo) ~ r" 1 Y^imi tytmt {k, r)Y tme (f) where 



nF e e ' v l n /'(E,E>)Y;, mti (k)exp 



i(kr- +Vi 9 ) 



(8) 



with ly? being the background phase shift in the ab- 
sence of interchannel and laser couplings. Now, 

equating the asymptotic form of r) with the 

boundary condition r -1 Y^e m [ s ^ n (kr — £tt/2 + 

m) - T eme (k)e l ( kr -<' 7 '/ 2+ ^}Y eme (r) we deduce the 

anisotropic scattering T-matrix element T( me (k) = 
E* m „ T imt ,e> mt , (E)Y* mei (k) where 



Te mt ,e'm e , (E) = T^ 9 S U 'S me m e , + ^Ft™* ( E > E ) e 



(9) 



with T\ g = —e^t" sin(r^ s ) being the T-matrix element in 
the absence of all laser and magnetic couplings. The 5-matrix 
element is obtained from Se me (k) = Ye me (k) — 2iTe me (k). 
The total scattering cross section can be obtained by means 
of the optical theorem as a tot = -^r Yltm e lm \ T hn t ,im e ], 
where the symbol Im stands for the imaginary part and g s = 
1(2) for two distinguishable (indistinguishable) atoms. The 
total elastic scattering cross section <j e i = J2e a & where at = 



Em, Vlmt and a £m t = "p 4 J2e'm e , \ T l 

inelastic scattering cross section erj„ e i = 
use of Eq. <j9j, we have 



m,/m ( ,| 2 - The total 
otot - del- Making 



4irg s 



E 

£rri£ 



+ E \< 



(10) 



For both the spontaneous linewidths 71 — > and 
72 — > 0, we have Ui ne i — > and the S ma- 
trix is unitary, that is m «,,r m( „ = 

The first term on the RHS of Eq.© is the isotropic part 
that does not depend on laser and magnetic fields. The sec- 
ond term proportional to Fj^ 1 ' (E, E) contains all the ef- 
fects of laser and magnetic fields on scattering. The term 

Etm/ (Ei E) given by Eq. © is a sum of three terms cor- 
responding to the three dressed bound-state amplitudes. For 
both lasers to influence p- or c?-wave scattering, it is essen- 
tial that the rotational quantum numbers of the two excited 
bound states should be such that both of them are allowed 
by selection rules for PA transitions from the p- or c?-wave 
ground scattering state. Once bound state 1 or 2 is populated 



10" 



10" 



ffl -12 

b 10 



10" 



1 
el 


r 

\ 




\ 

\ 

\ 




\ 

\ 


a. . 
inel 



50 



100 
9Ap 



150 



200 



FIG. 2: This figure is for 174 Yb atoms. Plotted are the total elastic 
scattering cross section tr e ( (solid) and total inelastic scattering cross 
section atnel (dashed) as a function of g2/"fs P for gi = 2007 sp at 
E = 100/iK. The value of ~f sp is 200 kHz. 



or their coherent superposition is formed, a nonzero partial- 
wave interaction in the ground continuum can be generated 
by two possible stimulated bound-free transition pathways. 
These pathways can interfere constructively under appropri- 
ate conditions resulting in enhanced higher partial-wave scat- 
tering amplitude. The sum of the two terms in F e-r ^ e ' (E, E) 
for £ ^ amounts to a coherent superposition of the ampli- 
tudes of those two transitions pathways. Furthermore, when a 
nonzero partial-wave ground scattering state is strongly cou- 
pled simultaneously to a pair of excited rovibrational states, a 
coherence between the two excited states due to the coupling 
is generated. In the next section, we discuss how this 
coherence influences ground-state collisional properties. 



III. RESULTS AND DISCUSSIONS 

The analysis presented in the preceding section is quite gen- 
eral since it is applicable to any ultracold atomic gas. Equa- 
tion ([3]) represents a nonperturbative eigenvalue solution of the 
model and is expressed as an anisotropic continuum dressed 
by two lasers and a magnetic filed. In the absence of the 
two lasers, the model reproduces the standard results of Fano- 
Feshbach resonances as can be verified from Eqs. ((9), (E), 
(O and (0. In the absence of one laser and the magnetic 
field, the model reduces to that of s-wave yD and (or) higher 
partial-wave OFR Jgt]. In what follows we apply this theoret- 
ical method to two specific atomic systems to show how the 
excited-state bound-bound coherence described by the term 
K-nn' affects ground-state elastic and inelastic scattering prop- 
erties. 

For numerical illustrations, we consider two atomic sys- 
tems: (A) fermionic 6 Li alkali atoms and (B) spin-singlet 
ground-state bosonic 174 Yb atoms. We assume that the two 
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PA lasers L\ and L2 are co-propagating and linearly polarized 
in the z axis for both systems. For this geometry, K nn i 7^ 
only for mi = Mi = M2 for both systems. 

A. Fermionic 6 Li atoms 

To select appropriate excited states of 6 Li2 for our nu- 
merical illustrations, we consider low-lying rotational levels 
of a relatively deeply bound vibrational state of electroni- 
cally excited 1 3 £+ potential which belongs to Hund's case b 
and asymptotically correlates to two separated (S1/2 + A/2) 
atoms. We calculate several ro-vibrational states supported by 
this potential. The purpose is to choose two excited bound 
states which lie far below the dissociation threshold of the ex- 
cited potential such that the possibility of laser transitions at 
the single-atom level can be safely ruled out. Furthermore, ro- 
tational energy spacings of such bound states should be large 
enough compared to other energy scales of our model. Hy- 
perfine interaction in such excited bound states is assumed to 
be negligible compared to the resonant dipole-dipole interac- 
tion between the two atoms. The laser couplings f2 x „ and 
^■isim depend on the respective Franck-Condon (FC) over- 
lap integral. To estimate FC integrals, we consider the v = 57 
vibrational state of the 1 3 S+ potential. This state has a vi- 
brational energy 2575.3 GHz below the dissociation threshold 
and an outer turning point around 30 do. The rotational energy 
spacings AE v j = E Vj j+i — E Vj j [where E Vt j stands for the 
binding energy of the (v, J) ro-vibrational state] are 3.4, 6.7 
and 10.4 GHz for J = 0, J = 1 and J = 2, respectively. The 
rotational quantum number is given by J = N = L + £ where 
L is the molecular electronic orbital angular momentum. The 
selection rule for laser transitions dictates that A A" = ±1. 

Exchange symmetry of a pair of spin polarized ground state 
fermionic 6 Li atoms allows only odd partial waves in the 
ground continuum. In such a case, p-wave scattering state of 
the ground continuum can be optically coupled to both J\ = 
and J2 = 2 rotational states by the two lasers. This can lead to 
the p-wave manipulation through optically generated coupling 
K nn i between the two excited rotational states. 

Let us now consider two-component fermionic 6 Li alkali 
atoms initially prepared in hyperfine spin f — 1/2 with two 
magnetic components mj = —m't = 1/2. In the presence of 
a strong magnetic field B (near B ~ B Q — 834. 1 G for MFR), 
neither atomic hyperfine quantum number / nor its projection 
nif along the quantization axis is a good quantum number, 
but the total magnetic quantum number Mp = mj + mp 
of the two atoms with hyperfine numbers / and /' remains a 
good quantum number at any magnetic field strength. In case 
of two-component 6 Li atoms with Mp = 0, dipole selection 
rules allow d-wave to be coupled to both Ji = 1 and J2 = 3 
rotational states by the two lasers while s-wave can be coupled 
to only J\ = 1 rotational state. Henceforth, we focus only on 
d-wave manipulation. 

To estimate ground-state input parameters, we calculate 
ground-state channel potentials in the presence of a magnetic 
field making use of the singlet and triplet ground-state poten- 
tial data available in the literature Il30l Bill . There are five 



asymptotic channels for M p = 0. We consider the lowest two 
asymptotic channels as a two-channel model for our numeri- 
cal work. The method of calculation of ground-state channel 
potentials is described in some detail in Appendix B. Since 
the closed channel quasibound state | \) nas s-wave rotational 
angular momentum, by selection rules we have fl xn ^ for 
n — 1 and f2 x „ = for n — 2. However, both the s-wave bare 
continuum and the | x) state can be indirectly coupled to the 
d-wave via the J = 1 bound state by two-photon process with 
both photons coming from laser L\ in the strong-coupling 
PA regime. Furthermore, the rotational selection rule dictates 
Agg* 2 = and hence K, nn ' = V„„< - iQ nn '/2. Since V nn > 
is a small shift between the two excited states, for our numer- 
ical work we set V nn ' = 0. Thus we have /C nn < = —iGnn' /2 
where Q nn , = 27rA^^A^™' 2fr ^' implying that this coupling 
term arises from the dipole interactions of the two excited 
bound states with the d-wave ground scattering state only. 

Let us suppose that the frequencies of both lasers are tuned 
such that E - M n - ££ hift ~ for n = 1, 2. Then under the 
conditions mentioned above, the the numerator on the RHS 
of Eq. (0]) becomes real provided the lasers do not introduce 
any phase as can be verified from Eqs. (1A251 and (IA26t . It 
then follows from Eq. ( [Tol l that the terms that make positive 
contributions to <7i ne i are proportional to {$$\ + ji)/\T>i\ 2 or 
(%+72)/|£> 2 | 2 , where V n = S n + *(Sf fl +g n )/2. Therefore, 
in order to suppress inelastic scattering, the terms Sf„ have to 
be reduced. Now, since <£ n = T™ f + Z ^>o,m „ r m me + Tw, 
and under the conditions r, m / of Eq.( |A36t becomes equal 
to — \Qnn'\ 2 /{Tn' + Jn') implying suppression of inelastic 
scattering due to optically induced coherence between the 
two excited bound states. Furthermore, since the d-wave 
(t = 2, mi) elastic cross section is proportional to the term 

X^' m , \Pi=2m t \ 2 w hi c h is basically the summation of am- 
plitudes for all the possible transitions from the two bound 
states to the d-wave ground continuum, the d-wave elastic 
scattering cross section would be enhanced when those transi- 
tion amplitudes add up coherently or in other words, interfere 
constructively. The numerical results discussed below further 
corroborate these analytical findings. 

For 6 Li atoms, the s-wave MFR is expected to influence 
the d-wave scattering process at relatively low energy. We find 
that for energies greater than 100 yuK, the MFR has practically 
no influence. Therefore, for higher energies d-wave manipu- 
lation occurs mainly due to two lasers. This means that for 
higher energies our method of d-wave manipulation is essen- 
tially all optical. This optical method can be applicable for a 
wide range of energies in the sub millikelvin regime. We ex- 
press the partial stimulated linewidth as V n ,tm t = Sn^n,(m,, 
where g n is the dipole coupling proportional to the intensity 
of the ?ith laser and J^f n ,im e is the angular part of the coupling. 
We assume that 71 = 72 = j S p and scale all energy quantities 
by hj sp . For all our numerical calculations, lasers Li and L2 
are assumed to be on resonance with PA transitions between 
the ground continuum and the excited bound states 1 and 2, re- 
spectively. For 6 Li, an energy-dependent s-wave background 
phase shift is estimated from the experimental parameters of 
the MFR reported in Ref. |4] as discussed in Appendix B. For 
6 Li atoms, the height of the d-wave centrifugal barrier is about 



6 



50 mK and the peak of the barrier is located at a separation of 
around 50ao. Therefore, the unperturbed d-wave scattering 
amplitude in the inner region of the barrier will be negligible 
unless the initial collision energy is relatively higher. 

In Fig.l, we show s-wave and d-wave elastic and total in- 
elastic scattering cross sections of 6 Li as a function of the 
intensity of the second laser for two different values of the 
intensity (gi) of the first laser at magnetic fields B near res- 
onant value Bq [ Figs. 1(a) and 1(c)] and far away from the 
MFR [Figs. 1(b) and 1(d)]. Comparing Fig. 1(a) with Fig. 
1(b) and Fig. 1(c) with Fig. 1(d), we infer that while the 
optical coupling of the MFR with the d-wave leads to sig- 
nificant enhancement of the d-wave scattering cross section 
in the strong-coupling regime, it has a marginal effect on in- 
elastic scattering. It can also be noticed that the intensity of 
the second laser has almost no effect on the s-wave scattering 
amplitude for a fixed value of B. Comparing Fig. 1(a) with 
Fig. 1(c), we observe that with the increase of the intensity of 
the first laser, the d-wave scattering cross section increases 
when B is tuned near the MFR. Since J2 = 3 rotational state 
is not accessible from the s-wave scattering state, the free- 
bound stimulated linewidth (r„) of the second (n = 2) bound 
state is almost entirely composed of d-wave partial stimu- 
lated linewidth. At 100 /iK temperature, in order to to make 
r2 ~ 1 MHz, the intensity of laser L2 has to be about 500 
kW/cm 2 , while to make T% ~ 1 MHz in microkelvin tem- 
perature regime the intensity of laser L\ is found to be of the 
order of 10 kW/cm 2 . The unperturbed or background d-wave 
phase shift at 100 /jK temperature is found to be smaller than 
that of the s-wave background phase shift by four orders of 
magnitude. This means that the unperturbed d-wave scatter- 
ing cross section in the microkelvin temperature regime is ex- 
tremely small. Since the centrifugal barrier height is inversely 
proportional to the reduced mass of the two atoms, for heavier 
atoms the d-wave barrier height is lower. Therefore, in case of 
heavier atoms it might be possible to manipulate the d-wave 
interaction more efficiently with much lower laser intensities 
as shown in the next section. 



B. Bosonic 174 Yb atoms 

The method of manipulation of interatomic interaction in 
174 Yb atomic gas is all optical because 174 Yb has no mag- 
netic moment. For d-wave manipulation, as in the case of two- 
component 6 Li atoms, it is necessary to tune the two lasers L\ 
and L2 near the exited rotational levels J\ = 1 and Jo = 3, 
respectively. For the excited Yb2 molecule, J = J e +£ where 
J e is the total electronic angular momentum. The two excited 
states are assumed to belong to 0+ (Hund's case c) molecu- 
lar symmetry, that is, the projection of J e on the internuclear 
axis is zero. The s- and d-wave background phase shifts and 
other input parameters for Yb are estimated from the data re- 
ported in Refs. ll28l l29tl . The ground-state potential of Yb 
is singlet only with no hyperfine structure and therefore the 
ground-state scattering of Yb is essentially single-channel. 

We plot total elastic and inelastic scattering cross sections 
of 174 Yb as a function of 92 in Fig. 2. From this figure, we 
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FIG. 3: These plots are for 174 Yb atoms, (a) ao (dashed) and 02 
(solid) are plotted as a function of g^/fsp- (b) <T2m 2 Vs. (?2/7s P for 
?7i2 = (solid), m,2 = 1 (dashed) and 7712 = 2 (dashed-dotted). The 
other parameters are the same as in Fig. 2. 



notice that for extreme strong-coupling regimes inelastic scat- 
tering can be suppressed by two orders of magnitude as com- 
pared to elastic scattering. 92 = 2OO7.jp corresponds to 40 
MHz. From the theoretical results reported in Ref. ll29ll . we 
estimate that at a laser intensity of 1 W/cm 2 , the d-wave stim- 
ulated linewidth of 174 Yb at 100 fxK energy is about 0. 1 MHz. 
Therefore, a 40-MHz stimulated linewidth will correspond to 
laser intensity of about 0.4 kW/cm 2 . Since the spacing be- 
tween the two excited bound states can be more than 50 MHz, 
the possibility of exciting multiple states by the laser 1 or laser 
2 is unlikely. 

We compare s- and d-wave scattering cross sections in Fig. 
3(a) which shows that the d-wave scattering cross section 
can exceed the s-wave one by an order of magnitude. Since 
fc>nn' 7^ on ly f° r m i — Mi = M2, tC nn i influences d-wave 
scattering for — 1 < m,(—2 < 1 only while me = ±2 remains 
unaffected by this term. Figure 3(b) demonstrates that <72 m2 
for 1712 = is higher compared to those for 7712 = 1 and 
?7i2 = 2 in the low 52 regime but <72m 2 f° r different 7712 values 
tend to merge and saturate in the high 92 regime. This means 
that in the low 32 regime the enhancement in d-wave scat- 
tering occurs predominantly with an interatomic axis oriented 
along z axis or the axis of quantization. 

From Fig. 2 we estimate the reduced inelastic rate of 174 Yb 
for a typical number density of 10 13 /cm 3 . Taking tji ne i ~ 
2 x 10 -13 cm 2 and the velocity at 100 /iK temperature to be 
~ 10 cm/s, we find an inelastic rate ~ 20 s~ 4 . The elastic rate 
for <r e ; ~ 10~ u cm 2 would be two orders of magnitude larger 
than the inelastic rate. These numbers indicate that optical 
manipulation of d-wave interaction in Yb is possible. 
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IV. CONCLUSIONS 

In conclusion, we have presented an analytical method for 
coherent manipulation of higher partial-wave atom-atom in- 
teraction by a pair of laser fields and an optional magnetic 
field in the strong photoassociative coupling regime. At low 
energy, this method is particularly applicable for manipula- 
tion of p- and (i-waves. Although, we have discussed p-wave 
manipulation, we have focused only on d-wave manipulation. 
The key effect predicted is the optically induced coherence be- 
tween two rotational excitations. As a proof-of-principle, we 
have demonstrated the significance of this coherence in sup- 
pression of inelastic scattering and enhancement of d-wave 
elastic scattering. We have derived an analytical expression 
for the ground-state scattering T matrix which clearly con- 
tains the effects of quantum interference. An analytical form 
of the two-body T matrix is important as it may serve as a key 
ingredient for any theoretical treatment of interacting many- 
particle systems. In this work, we have discussed application 
of our method to two-component fermionic 6 Li and bosonic 
17 Yb atoms. However, our method is applicable to all the 
atomic species that are currently being used in cold atom 
research. It would be particularly interesting to apply our 
method to explore manipulation of p- or d-wave interaction 
in heavier fermionic atoms such as 171 Yb, 173 Yb and 87 Sr, all 
of which have very narrow linewidths. In order to make this 
method widely applicable, it is also important to find ways of 
reducing the natural linewidth of molecules by either tailor- 
ing vacuum fluctuations or by an interplay between light- and 
vacuum-induced coherences ll32ll . 



Schrodinger equation H \ E, k) = E \ E, k), making 

use of the expansions Be = Ylt'm j ^E m ''^*m ( ,(^) an d 

A nE = Y,e>m e , A e nE e 'Y e *, m ^ (*0> we obtain the following set 
of coupled algebraic equations 

x <fe(i?)> (Al) 

(E x - E)B^' + £ Ml'?' = - / dE>V x0 (E>) 

x c£%, (A2) 

and 

(E'-E)C^(E) + £A&f 

n,M„ 

+ V 0x {E')B E mi '5 m = Q (A3) 

where S n = E n /H — (ljl„ — ^>a) with E n being the binding 
energy of n-th excited bound state measured from the thresh- 
old of the excited state potential, Wi n is the laser frequency 
of n-th laser and oja the atomic transition frequency. Since 

£' m 

the coefficient C E , e ' n( is required to fulfill scattering bound- 
ary conditions at large separation of the two atoms, we can 
express Eq. iA3i in the form 



Acknowledgments n i'm e , (J? ^ _ xfTP j?i\x x ^Ox(E') u e'm 



Ce'CSP) = 5(E-E>)6 i , i ,5e mi , + j^B E m *'5 l 



The author is thankful to C. Salomon, R. Hulet and Saikat _ Ai n Mn (E') ,, 

y im t \ ) A tm t , (A4) 

Z_^i E — E' 

n,M„ 



Ghosh for discussions and comments. 



Appendix A: DERIVATION 

Here we present detailed derivation of the anisotropic 
dressed continuum of Eq. (|3). From the time-independent 



The partial-wave symbols appearing in the superscript and 

subscript of C E ™^ 1( (E) refer to the incident and scattered par- 
tial waves, respectively. Substituting Eq. (IA41 > in Eqs. ( lAlb 
and (1A2K we obtain 



(s n - e)4™" + n nx B l E m * = -aJX(jso- 



dE 



,V 0x (E')^ Mn (E') 
E-E 1 



B 



t' m^t 



E 

£rri£ 



dE 



JA ,tf"(^)| 2 
E-E' 



x a: 



EE 



M^HE')A^ M JE>) 



dE 



E-E' 



A e'm e , 



(A5) 



(E x - E)B i E m " + £ n xn Ai7' = -V x0 (E)5 e , - [dE 



I \ Vxo( E ')\ 2 B e'me> 

E-E' E 



Y^jdE 

l, Mr, 



, V x0 (E')A S e ' mtl 

E-E' Ke (A6) 



r 



where n! ^ n. In writing the third term on the right side imt)(lmi \ e ■ d \ JM') = ^2 M i(JM e ■ d \ £mi)(£me \ 



e • d | JM')6 M ,M' = \{^fnt [ e - rf | JM)| 2 , where e stands 
for laser polarization and a! is a unit vector pointing toward the 
molecular dipole moment. Taking E — > E + ir\ with 77 = + , 
we have 



^^OA^O = vt o )[E) _ l7:Vox{E) 

x A^ Mn (^) (A7) 
where 



E-E' 



is an effective interaction between the closed-channel bound 
state I x) an d me excited bound state | b n ) mediated through 
the s-wave part of the ground continuum. Similarly, we have 



dE > \ V 0x(E')\ 2 =E *m_ i £f 
E-E' x 2 ' 



(A9) 



where 



E shift = v I ,\Vo x (E')\ 2 
x I E-E 1 



(A12) 



is the shift of the closed-channel bound state for its interac- 
tion with s-wave part of the open channel continuum, Tf = 
2tt\V[) X {E)\ 2 is the magnetic Feshbach resonance line width, 

«f. — 'P j dE' — 



E _ E , is the partial light shift and 

F n £ me {E) = 2n\Ag^f n (E)\ 2 is the partial stimulated line 
width of the n-th excited bound state. Here 



and 



jw|A^ (E )l _ pshift .r,i£m f 



(A10) 



E-E 1 



Q^ ] =2,K^{E)KTJmJE) 



JnM„ ( 771 \ a <m< 



(A14) 



Substituting Eqs.dATb, dA9l ), dAlOb and dAl lb into Eqs. 



(All) (|A5T > and (lAol we obtain 



(£-*n-2i^ + *T n /2)4£* -{^^ (A15) 



(2? - E x - Ef* + iT f /2)B e E m *' i Q xn + K ( f 0) - i^ M -V x0 {E)}K7' = V x o(E)S m (A16) 

I 



.,M„ 



where E^ = J2em e ^Vem f i s t ' le tota l light shift due to pressed in a compact form 
laser L n , T n = J2i me ^ntmt i s me corresponding total stim- 
ulated linewidth and 



B 



^=Ete f) -4^) (A17) 



(e + i)r/ 
x 7rA^ M «y x0 (^)4>'] 



F x0 (£)^<o+ E (9n/-0 

(A18) 



.A/„ 



where 



E-E x - E s x hilt 
T772 



(A19) 



is an effective complex coupling between the two excited 
bound states | b n ) and | b n >^ n ) induced by the two lasers 

through optical dipole interactions. The superscript LL is used is the dimensionless collision energy measured from the 
to indicate the use of two lasers. Equation dAl 6b can be ex- shifted binding energy E x + E^ llt of the closed channel qua- 
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sibound state | \) ar, d 

q n f = 



il + V {1=0) 
nA J oS M »V xQ (E) 



(A20) 



is the Fano-Feshbach asymmetry parameter that describes 
quantum interference between two transition pathways (i) 



L„ La 



x) an d (ii) I b n ) Ln LahC1 > E£ — 0m ( 



I K) 

0)bare — >\ X}- Substituting Eq. (|A18l l into Eq. ( IA15I ), we 
have 



E-S n 



(q n f -i) 2 T n 



no 



A* 



q n f - i 

e + i 



I'm,, 



A- Mn (S)^ + A^X(S) 



Here 



where 



JC nn ' (M„, M n >) = JC^, + /C 



LL 

an/ 



, J„,M„ 



(A21) 



(A22) 



(A23) 



Note that, when both the lasers are linearly polarized along 
the quantization axis, we have K. n n'(M n ,M n ') = unless 
M„ = M n i. Putting n = 1, 2 and n' = 1, 2 with n' ^ n in 
Eq.(A17), one can write down two coupled equations which 
can be solved to obtain 



A t'm t i I „ f \ _ MlO^'.O + (1 — <^'o) 
^nB l iV2 nJ - p 



where 



(A24) 



(A25) 



(A26) 



and 



00 



e + i 



(A27) 



In Eq. dA251 >. /3„ e is a dimensionless factor defined by 

a _ ?»/ + e 



(A28) 



e + z 

and in Eq. ( |A27t (E) is given by 

=E-hS n - ££ hift + ihfr n + T n {E)}/2 (A29) 

In the above equations n' ^ n and M n = M n i . We have here 
phenomenologically introduced the spontaneous linewidth j n 



by making the replacement T„ — > T n +j n . Writing Eq. (|A271 > 
in the form 



V n (E) = g n + ih(& n + 7„)/2, 



we have 



$ n — E — hS n — £^ hlft — ff ? " lf — S nn i 



pmf 
1 nOO 



e>o,m e 



(A30) 



(A31) 



(A32) 



where 



mf 



(q 2 



nf 



l)e - 2g n/ 



nOO 



^n00 — l/3rie| 2 r 



nOO 



2 

r?m' — ^nn'^n'n\ 



(A33) 

(A34) 
(A35) 

(A36) 



Appendix B: GROUND STATE OF TWO-COMPONENT b Li 
ATOMS IN A MAGNETIC FIELD 



In the presence of a magnetic field, the angular basis states 
of the form | //', m/, mf, Mp) (diabatic basis) of two sep- 
arated atoms for different values of /, /', m/ and ntf get 
mixed up. The Hamiltonian of two ground state atoms at large 
separation in the presence of an external magnetic field is then 
diagonalizable in these diabatic basis and the eigenfunctions 
are regarded as the scattering channels. However, these are not 
diagonal at short separations where the Hamiltonian is diag- 
onalizable in the molecular adiabatic basis. In the absence of 
hyperfine interactions, there are two adiabatic molecular po- 
tentials corresponding to the spin singlet and triplet, both of 
which asymptotically go as ~ — 1/r 6 . The diabatic and adi- 
abatic basis are related via transformation involving Wigner 
9-j symbol. For Mp = there are five channels which be- 
come coupled at intermediate separations ( < 30ao, where ao 
is Bohr radius). The coupling is most significant at separa- 
tions around 20 do. We consider the two lowest asymptotic 
channels with Mp = 0. The lower channel | g{) is open 
and the upper channel | 32) is closed. The asymptotic en- 
ergy spacing between these two channels at Bq is calculated 
to be 2.25 GHz. The open channel is mostly electronic spin 
triplet (about 90% triplet) while the closed channel is largely 
spin singlet. Exchange symmetry of a pair of two-component 
fermionic atoms with Mp = allows only even partial waves 
of relative motion of the pair, while for a pair of single compo- 
nent fermionic atoms only odd partial waves are allowed. For 
6 Li atomic gas initially prepared in a 50:50 mixture of two 
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hyperfine components of / = 1/2 and Mp = 0, even partial- 
wave atom-atom interactions are expected to be predominant. 

At low energy, the Feshbach resonance width T/ is given 
by Tf/2 ~ kcibgTo, where k is the wave number re- 
lated to the collision energy E = H 2 k 2 /(2m) with m be- 
ing the reduced mass of the two atoms, a,b g = — 1405ao 
is the back ground scattering length, Tq is a parameter re- 
lated to the width A = —300 G of the zero crossing by 
A = J 2 -, with 5fi = 2/ib (where [Lq is the Bohr mag- 
neton) being the difference of magnetic moment of closed- 
channel quasibound state from that of the two free atoms 
in the open channel. The energy-dependence of the Fesh- 
bach resonance phase shift rj res is given by — cot ij res = 
e = E f Ei ^ - E ffz h , where E th is the threshold of the 



open channel that depends on the applied magnetic field, and 
E c — E t h = 5(1(13 — Bq). The background s-wave phase 
shift (with effective range of background open channel poten- 
tial assumed to be zero) t]q 9 is given by cot^ 9 ~ — j^-- 
Sine the coupling between open and closed channels occurs at 
an intermediate separation (kr 1), we can write Tf/2 — 
7T\V x0 \ 2 = n\fdri, x (r)V(r)4, E j= (r)\ 2 ~ IT^F 1 ^ 
where V(r) stands for coupling between the two channels 
and ipE,e=o(r) = (r | E,£ = 0) bare is the energy nor- 
malized background (open channel) scattering state, which 
asymptotically behaves as ipE,oo(r — > oo) ~ sin(£;r + r^ 9 ) 
and ip x {r) = (r\x) 
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